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Abstract 

The /5t-differential inclusive production cross sections of the prompt charmed mesons D", D+, and 
D*+ in the rapidity range \y\ < 0.5 were measured in proton-proton collisions at ^ = 7 TeV at 
the LHC using the ALICE detector Reconstructing the decays D"— >K";r+, D+ K"7r+;r+, 
D*+ -> 0"^:+, and their charge conjugates, about 8,400 D", 2,900 D+, and 2,600 D*+ mesons with 
1 < Pt < 24 GeV/c were counted, after selection cuts, in a data sample of 3.14 x 10^ events collected 
with a minimum-bias trigger (integrated luminosity Lint = 5 nb^'). The results are described within 
uncertainties by predictions based on perturbative QCD. 



*See Appendix lAlfor the list of collaboration members 
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1 Introduction 

The study of the production of hadrons containing heavy quarks, i.e. charm and beauty, in proton-proton 
(pp) collisions at LHC energies provides a way to test, in a new energy domain, calculations of quantum 
chromodynamics (QCD) processes based on the factorization approach. In this scheme, cross sections are 
computed as a convolution of three terms: the parton distribution functions of the incoming protons, the 
partonic hard scattering cross section calculated as a perturbative series in the strong interaction coupling 
constant, and the fragmentation function, parametrizing the relative production yield and momentum 
distribution for a charm quark hadronizing to a particular species of D meson. Recent implementations 
of such calculations, at the perturbation level of next-to-leading order or at fixed order with next-to- 
leading-log resummation (FONLL) HI describe well the beauty production cross section measured in 
pp collisions at y/s = 1.96 TeV at the FNAL Tevatron collider |l2]-@l and in pp collisions at y(s = 7 TeV 
at the CERN Large Hadron Collider (LHC) ||5]|6l. The production cross section of charmed hadrons (D 
mesons) at Tevatron is reproduced within the theoretical uncertainties of the calculations as well ||7]^. 
However, the comparison suggests that charm production is underestimated by the results obtained with 
the central values of the calculation parameters, as observed also in pp collisions at the BNL RHIC 
collider at the lower energy of y/s = 200 GeV ifTOlfTTl . In this context, it is particularly interesting to 
perform the comparison for charm production at the LHC energy, which is more than three times higher 
than at the Tevatron. Furthermore, at LHC energies, the measurement of charm production in the low 
transverse momentum {pi) region probes the parton distribution functions of the proton at small values 
of parton fractional momentum x and squared momentum transfer Q^. For illustration, using a simplified 

2 — )• 2 kinematics at leading order, c quarks {m^ w 1.5 GeV/c^) with pi^2 GeV/c and rapidity y ~ 
probe the parton distribution functions at ~ 7 x 10^^ and ~ (5 GeV)^, where the gluon component 
is dominant. In this kinematic regime, the gluon distribution may reach the level of saturation, leading to 
a measurable departure of the observed cross sections from the expectations based on the factorization 
approach (see e.g. fV2\ ). 

We report on the measurement of the production cross section of the prompt charmed mesons D", D+, 
and D*+, in pp collisions at ^/s = 7 TeV, reconstructed in the range \ < pi <2A (16 for D°) GeV/c and 
\y\ < 0.5 with the ALICE experimental apparatus |[T3]| . using data collected in 2010. The detector layout 
and the data sample used for the measurement are described in section |2] The D meson reconstruction 
procedure, the raw yield extraction, and the corrections applied to obtain the production cross sections, 
including the estimation of the prompt D meson fraction, are presented in sections [3] and ID Finally, the 
resulting pt-differential cross sections are reported in section 1431 and compared to QCD predictions. 

2 Detector layout and data collection 

D", D+, and D*+ mesons, and their charge conjugates, are reconstructed from their decays into charged 
hadrons in the central rapidity region utilizing the tracking and particle identification capabilities of the 
ALICE central baiTcl detectors. 

The ALICE apparatus is described in detail in |[T3l . It consists of a central barrel covering the pseu- 
dorapidity interval < 0.9, a forward muon spectrometer, and a set of small detectors in the forwai^d 
regions for trigger and event characterization purposes. Only the detector features that are relevant for 
the D meson analysis are presented here. The ALICE global coordinate system is right-handed, with the 
origin coinciding with the geometrical centre of the central barrel, the z axis directed along the beam 
line, the x axis in the LHC (horizontal) plane, pointing towards the centre of the accelerator, and the y 
axis pointing upward. The central barrel detectors are positioned within a large solenoid magnet, with a 
field B = 0.5 T, pai^allel to the beam line. Data collected with both magnet polarities were used for this 
analysis. 

The innermost detector of the ALICE central barrel is the Inner Tracking System (ITS) which is made 
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of six cylindrical layers of silicon detectors, with radial distance to the beam-line between 3.9 cm and 
43.0 cm. The two innermost layers, with average radii of 3.9 cm (about 1 cm from the beam vacuum 
tube) and 7.6 cm, are equipped with Silicon Pixel Detectors (SPD), comprising 9.8x10^ pixels of size 
50 (r^) X 425 (z) /xm^, with intrinsic spatial resolution of 12 (r^) x 100 (z) /im^. The signals of the 1,200 
SPD readout chips provide a fast trigger through a programmable logic. The two intermediate layers, 
at radii of 15.0 and 23.9 cm, are made of Silicon Drift Detectors (SDD). They allow one to measure 
the hit position along z with resolution better than 30 /im from the centroid of signals collected on the 
anodes, and to determine the r0 coordinate from the drift time with a resolution that depends on the 
level of calibration, as discussed below. The two outermost layers are made of Silicon Strip Detectors 
(SSD), located at radii of 38.0 and 43.0 cm, consisting of double-sided silicon strip sensor modules, 
with an intrinsic spatial resolution of 20 (r^) x 830 (z) ^um^. The total material budget of the ITS is 
on average 1.1% of radiation length for tracks crossing the ITS perpendicularly to the detector surfaces 
(t] = 0) II13II14II . These features enable measurement of the track impact parameter (i.e. the distance of 
closest approach of the track to the primary interaction vertex) in the bending plane (r0) with a resolution 
better than 75 pLm for transverse momenta Pt > 1 GeV/c (see section [TTI ). 

The ITS modules were aligned using survey information, cosmic-ray tracks, and pp data, with the 
methods described in |[T4l . For the residual misalignment along the r0 coordinate, a r.m.s. of about 
8 /im for SPD and 15 pLra for SSD modules was estimated II14|]151 . For SDD, with the current calibration 
level, the space point resolution along r0 is 60 /im for those modules (about 60% of the total) that 
do not suffer from significant drift field non-uniformities. To account for the fact that a correction for 
these effects was not applied in the reconstruction of 2010 data, a systematic uncertainty of 300 /xm 
was assumed for SDD points. Along z, the estimated precision of the alignment is 50 /im for SPD and 
SDD and a few hundred jUm for SSD. These values are added in quadrature to the uncertainty on the 
reconstructed ITS hits in the track reconstruction algorithm. In the detector simulation, to account in 
an effective way for the residual misalignment, the ITS modules ai^e randomly displaced with respect to 
their ideal positions according to the estimated precision of the alignment. 

At larger radii (85 <r < 247 cm), the cylindrical Time Projection Chamber (TPC) |[T6l provides track 
reconstruction with up to 159 three-dimensional space points per track, as well as particle identification 
via the measurement of the specific energy deposit AE/dx. The TPC has an active length of 500 cm 
along the z direction and its 90 m^ gas volume is filled with a mixture composed of Ne (85.7%), CO2 
(9.5%), and N2 (4.8%). Its position resolution is 1100-1250 jira along the z axis (corresponding to the 
drift direction) and 800-1 100 pLm along r0. Using cosmic-ray muons and data taken in pp collisions, the 
relative dE'/dx resolution was measured to be about 5.5% for tracks that cross the entire detector flQ . 

The charged particle identification capability of the TPC is supplemented by the Time-Of-Flight detector 
(TOF), that is based on Multi-gap Resistive Plate Chambers (MRPCs) in a cylindrical configuration at 
radius 370-399 cm from the beam axis, with readout consisting of 152,928 sensitive pads of dimension 
2.5 X 3.5 cm^. The TOF resolution on the particle arrival time is at present better than 100 ps. The 
start time of the collision (event time zero) is measured by the TO detector, an array of Cherenkov 
counters located at +350 cm and —70 cm along the beam-line, or, for the events in which the TO signal 
is not present, estimated using the particle arrival times at the TOF. The particle identification is based 
on the difference between the measured time-of-flight and its expected value, computed for each mass 
hypothesis from the track momentum and length. The overall resolution on this difference is about 
160 ps. In this analysis, the time-of-flight measurement provides kaon/pion separation up to a momentum 
of about 1 .5 GeV/c. Results from the TOF commissioning with cosmic-ray particles are reported in ifTTl . 

The data sample used for the analysis presented here consists of 314 million minimum-bias (MB) 
events, corresponding to an integrated luminosity Ljnt = 5 nb^\ collected during the 2010 LHC run 
with pp collisions at y/s = 1 TeV. The minimum-bias trigger was based on the SPD and VZERO 
detectors. The VZERO detector consists of two arrays of 32 scintillators each, placed around the beam 
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vacuum tube on either side of the interaction region. The two aiTays cover the pseudorapidity ranges 
—3.7 < T] < — 1.7 and 2.8 < T] < 5.1. Minimum-bias collisions were triggered by requiring at least 
one hit in either of the VZERO counters or in the SPD (|t7| < 2), in coincidence with the arrival of 
proton bunches from both directions. This trigger was estimated to be sensitive to about 87% of the pp 
inelastic cross section lITSl . It was verified on Monte Carlo simulations based on the PYTHIA 6.4.21 
event generator |[T9l (with Perugia-0 tuning |[20l ) that the minimum-bias trigger is 100% efficient for D 
mesons with pt> I GeV/c and \y\ < 0.5. Contamination from beam-induced background was rejected 
offline using the timing information from the VZERO and the con^elation between the number of hits and 
track segments (tracklets) in the SPD detector. The instantaneous luminosity in the ALICE experiment 
was limited to 0.6-1.2 x 10^^ cm^^s^^ by displacing the beams in the transverse plane by 3.8 times the 
r.m.s of their transverse profile. In this way, the interaction probability per bunch crossing was kept in the 
range 0.04-0.08, with probability of collision pile-up below 4% per triggered event. The luminous region 
was measured with high precision from the distribution of the interaction vertices reconstructed from the 
charged particles tracked in the ALICE central baiTcl detectors, yielding a'"™'"°"** « (jiummous 
50 jxm in the transverse plane and (7_i™i"°"s ~ 4_6 cm along the beam direction (the quoted ranges are 
due to the variations of the beam conditions during the data taking). Only events with interaction vertex 
in the range |z| < 10 cm were used for the analysis. 



3 Reconstruction of D meson decays 

The measurement of charm production was performed by reconstructing three decay modes for D 
mesons, D° K-7r+ (with branching ratio, BR, of 3.87 ±0.05%), D+ ^-71+11+ (BR of 9.13 ± 
0.19%), and D*(2010)+ D°7r+ (strong decay with BR of 67.7 ±0.5%) with D° K"7r+, together 
with their charge conjugates II211 . The D*' and D+ mesons have mean proper decay lengths cT 123 
and 312 pLra, respectively ll2ll . Their decay secondary vertices are therefore typically displaced by a 
few hundred /im from the primary vertex of the pp interaction. The analysis strategy for the extraction 
of the D" and D+ signals from the large combinatorial background due to uncorrelated tracks is based 
on the reconstruction and selection of secondary vertex topologies that have significant separation from 
the primary vertex. The identification of the charged kaons in the TPC and TOE detectors provides 
additional background rejection in the low-momentum region. A particle identification strategy that has 
high efficiency for the D meson signal (see section |4~TI ) was adopted. Einally, an invariant mass analysis 
was used to extract the signal yield. In the D*+ — D°;r+ case, the decay vertex cannot be resolved from 
the primary vertex. The analysis exploits topological selections on the D°, together with the sharp peak 
in the difference between the invariant mass of the three final state hadrons and that of the two D° decay 
prongs. Given that the mass difference Am = 171^,*+ — m^o 145.4 MeV/c^ 11211 is only slightly larger 
than the charged pion mass, for low pi D*+ mesons the produced pion has low momentum and is referred 
to here as a 'soft pion'. 

3.1 Track and vertex reconstruction 

The procedure for track reconstruction in the central barrel detectors II13II22I1 starts with a first determina- 
tion of the primary vertex position, performed by correlating hit pairs (tracklets) in the two layers of the 
Silicon Pixel Detector. The same algorithm is used to search for vertices from pile-up collisions starting 
from the tracklets that do not point to the first found vertex. An event is rejected from the analyzed data 
sample if a second interaction vertex is found, it has at least 3 associated tracklets, and it is separated 
from the first one by more than 8 mm. The remaining undetected pile-up is negUgible for the analysis 
described in this paper. Track seeds are built using this first estimate of the interaction vertex position 
together with pairs of reconstructed space points in adjacent readout pad rows of the TPC. Tracks are 
then projected inward in the radial direction using a Kalman filter algorithm, incorporating space points 
in the TPC and then hits in the six layers of the ITS (referred to here as TPC-i-ITS tracks). The ITS hits 
not associated to TPC-i-ITS tracks by this procedure are then used to search for ITS -only tracks, most 
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Figure 1: Track impact parameter (do) resolution in the transverse plane {r(j) direction) as a function of pt 
comparing data and simulation. This resolution includes the uncertainty in the primary vertex position, which 
is reconstructed excluding the track being probed. 



importantly to recover pions with SO < pt < 200 MeV/c that have very low reconstruction efficiency in 
the TPC, for geometrical reasons. ITS-only hacks are found using a hit grouping algorithm that projects 
angular windows radially outward from the primary vertex. The TPC+ITS tracks are then propagated 
outward in order to associate signals in the large-radius detectors that perform particle identification. Fi- 
nally, all tracks are re-propagated with the Kalman filter in the inward direction. The relative pi resolution 
at the primary vertex for this procedure is about 1% (6%) at 1 GeV/c for TPC-i-ITS (ITS-only) tracks. 
The last step of event reconstruction is the re-determination of the primary vertex position from the ac- 
cepted tracks The primary vertex coordinates and covariance matrix are obtained via an analytic 
minimization method applied to the tracks approximated by straight lines after propagation to their 
common point of closest approach. The algorithm is then repeated excluding the tracks with distance to 
the primary vertex, normalized to its estimated uncertainty, larger than 3, which are incompatible with 
being produced by primary particles. The vertex fit is constrained in the transverse plane using the in- 
formation on the position and spread of the luminous region. The latter is determined from the 
distribution of primary vertices averaged over the run and is tabulated as a function of time during the full 
data-taking period. The position resolution of the primary vertex reconstructed from tracks depends on 
the particle multiplicity. It was measured to be (Mm) « ^^^/^txackiets longitudinal direction and 
CJv.y (/im) ^ niin(a'"™"°"%600/A'^j°^^ygjJ in the transverse coordinates by fitting its dependence on the 
number of tracklets in the SPD (A'^nackiets). which con^esponds to about twice the multiplicity of primary 
charged particles per unit of rapidity. Thus, for the average luminous region spread al"™"™** w 40 ^m, 
the transverse position of the vertex has a resolution that ranges from 40 /^m in low-multiplicity events 
(i.e. below 10 charged particles per unit of rapidity) to about 10 /^m in events with a multiplicity of about 
40. 

Secondary vertices of D° and D+ meson candidates were reconstructed using tracks having |t]| < 0.8, 
Pt > 0.4 GeV/c, at least 70 associated space points (out of a maximum of 159) and ;^^/ndf < 2 in 
the TPC, and at least one hit in either of the two layers of the SPD. For tracks passing this selection, 
the average number of hits in the six ITS layers is 4.5-4.7, depending on the data taking period. This 
quantity is influenced by the fraction of inactive channels and its distribution is well reproduced in Monte 
Carlo simulations. For the D*+ soft pion, all TPC-i-ITS and ITS-only tracks with at least 4 hits in the ITS, 
including at least one in the SPD, and pi > 80 MeV/c were considered. Figure [T] shows the transverse 
momentum dependence of the transverse (r0) impact parameter (do) resolution achieved with the present 
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ITS alignment precision for tracks that satisfy the TPC and ITS selection criteria, for data and Monte 
Carlo simulations. The simulations utiUze GEANT3 |[24l and incorporate a detailed description of the 
detector material, geometry and response. Proton-proton collisions were simulated using the PYTHIA 
6.4.21 event generator llT9l with Perugia-0 tuning |[20i . The impact parameter resolution was estimated 
by fitting the inclusive distribution of do with respect to the event primary vertex, in intervals of pi. 
The fit function is the sum of a Gaussian, that accounts for the component due to prompt particles 
produced at the primary vertex, and two exponential functions, that account for secondary particles, 
mainly from weak decays of strange hadrons. The width a of the Gaussian provides an estimate of 
the do resolution, which includes the resolution of the track parameters and the primary vertex position. 
In order to obtain an unbiased estimate of d^, the primary vertex is recalculated excluding the track 
being probed. The figure shows that the do resolution measured in data, with values of 75 (20) pLra 
at pt = 1(15) GeV/c, is reproduced within about 10% by the Monte Carlo simulation incorporating 
the residual ITS misalignment described in section |2] The effect of the difference between data and 
simulation on the results of the D meson analysis is discussed in section l42l 

3.2 D meson selection 

D'', D+, and D*+ candidates were filtered by applying kinematical and topological cuts, and particle 
identification criteria. A fiducial acceptance cut IjdI <3'fid(A) was applied, with Vfid smoothly increasing 
from 0.5 to 0.8 in < pt < 5 GeV/c and = 0.8 above 5 GeV/c. For D" and D+ decays, the secondary 
vertex was reconstructed with the same algorithm used to compute the primary vertex from tracks. The 
resolution on the position of D° and D+ decay vertices was estimated by Monte Cai^lo simulations to be 
of the order of 100 pLm with little pt-dependence for pt > 1 GeV/c [23J. For the D*' and D+ selection, 
the primary vertex was recalculated for each D candidate, excluding the decay tracks. 

The cut variables for the three mesons are described in the following. The actual cut values are pt 
dependent and were tuned to optimize the statistical significance of the signal, resulting in a selection 
efficiency that increases with increasing pt. The cut values applied for D mesons at low pt are reported 
for reference in the next paragraphs. 

For D° mesons, the two decay candidate tracks were further selected with pt > 0.7 GeV/c (pt > 
0.4 GeV/c for 1 < pP < 2 GeV/c) and r0 impact parameter significance \dQ\/OdQ > 0.5. Secondary 
vertices were required to have a minimum displacement of 100 /^m from the primary vertex and a 
maximum distance of closest approach between the two tracks of 300 /^m. The cut | cos d*\ < 0.8, where 
0* is the angle between the kaon momentum in the D" rest frame and the boost direction, was applied 
to reduce the contamination of background candidates that do not represent real two-body decays and 
typically have large values of | cos 6*\. Well-displaced — K^n^ topologies are characterized by large 
and opposite-sign values of the decay track r(p impact parameters (d^ and d^) and good pointing of the 
reconstructed D° momentum to the primary vertex, i.e. a small value of the pointing angle ^pointing 
between the momentum and flight line. Due to the strong correlation of these two features in the signal, 
the requirement d^ xd^ < —{120 /^m)^ and cosQpomting > 0.8 was found to be effective in increasing 
the signal-to-background ratio. 

The D+ selection is based on a similar strategy to that for D° mesons. A looser cut on the pt of decay 
tracks, 0.4 GeV/c, was applied, due to the lower average momentum of the products of a three-body 
decay. The candidate triplets were selected on the basis of the sum of the distances of the decay tracks 
to the reconstructed decay vertex, the decay length, the cosine of the pointing angle, and the sum of the 
squares of the impact parameters of the three tracks. Typical cut values for lov/-pi candidates are: 
decay length larger than 800 /^m, cos ^pointing > 0.95, and sum of the squares of the three decay tracks 
impact parameters ZJq > (750 /xm)^. The topological selection cuts are, in general, tighter than for 
the D° case due to the larger cT of the D+ meson, resulting in a better separation between primary and 
secondary vertices, and to the higher combinatorial background in the three-particle final state. 
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Figure 2: Invariant mass distributions for D*^ (top) and D+ (middle) candidates, and mass difference distribution 
for D*+ candidates (bottom), for three p•^ intervals. The curves show the fit functions as described in the text. The 
values of mean (/i) and width (a) of the signal peak are reported (for D" and D+ they are expressed in GeV/c^). 



The D*+ candidates were filtered by applying kinematical selections on the final decay products and cuts 
on the topology of the D° decay. The single track minimum transverse momentum was set to 0.4 GeV/c 
for the D° decay tracks and 80 MeV / c for the soft pion track. The variables used to select the topology 
of the D'' decay are the same as for the D° analysis described above. However, a selection with higher 
efficiency could be applied in this case, because the background in the region around Am 145 MeV/c^, 
which is close to the phase space boundary, is much lower than that around the D° mass. In particular, 
for D*+ candidates with pt> 6 GeV/c, the topological cuts could be opened so as to select about 90% 
of the signal passing single track cuts. 

The particle identification selection used the specific energy deposit and the time-of-flight from the 
TPC and TOP detectors, respectively. In order to assign the kaon or pion mass to the decay tracks, 
compatibility cuts were applied to the difference between the measured and expected signals. For both 
dE/dx and time-of-flight, a 3 a compatibility cut was used. Tracks without a TOF signal were identified 
using only the TPC information, and tracks with incompatible TOF and TPC indications were treated as 
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Table 1: Measured raw yields for D , D+, and D*+ mesons, and their anti-particles, in a minimum-bias pp sample 
corresponding to 5 nb^' at ^/s — 7 TeV, in transverse momentum intervals. The systematic uncertainty estimation 
is described in section l4~2l 



Pi interval 


DO + d" 


A'' ± Stat. ± syst. 




(GeV/c) 




D+±D- 




_|_ Q*- 


1-2 


1531 ±233 ±340 


122±23± 


30 


78±16± 8 


2-3 


1978 ±168 ±190 


390±57± 


97 


244±26±10 


3-4 


1950 ± 129 ± 


75 


405 ±40 ±101 


363±29±11 


4-5 


1184± 78± 


40 


516±38± 


46 


478 ±36 ±14 


5-6 


623 ± 50 ± 


25 


361±31± 


33 


374±28±18 


6-7 


339 ± 32 ± 


13 


294±30± 


15 


279 ± 22 ± 14 


7-8 


199 ± 25 ± 


14 


213±27± 


11 


170±19± 9 


8-12 


427 ± 38 ± 


30 


434±30± 


22 


408 ±28 ±18 


12-16 


139 ± 27 ± 


14 


89±20± 


9 


115±12±10 


16-24 






52±14± 


5 


41± 6± 8 



non-identified, but still used in the analysis and considered to be compatible with both a pion and a kaon. 

n — 

Two-prong candidates were accepted (as D", D , or both) or rejected, according to the compatibility with 
the K^TT* final state. For D" candidates used in D*+ reconstruction, compatibility with the appropriate 
final state, consistent with the soft pion charge, was required. Particle identification was not applied to 
the soft pion tracks. In the case of the D+ — K^;r+7r+ decay, the particle with the opposite charge sign 
with respect to the D meson is a kaon. Hence, the triplets were rejected if the opposite sign track was 
not compatible with the kaon hypothesis, or at least one of the two same sign tracks was not compatible 
with the pion hypothesis. For all three D meson species, a comparison of the invariant mass distributions 
obtained without and with particle identification shows that this selection reduces the combinatorial 
background by a factor 2-3 in the low pt region, while preserving close to 100% of the D meson signal. 

The raw signal yields were extracted in the intervals of px listed in Table [TJ by a fit to the invariant mass 
distributions (or mass difference for the D*+), as shown in Fig. [2]for three selected pi intervals. For D° 
and D+ mesons, the fitting function consists of a Gaussian describing the signal and an exponential term 
for the background. In the D*' case, the contribution of signal candidates with wrong mass assignment 
to the final state hadrons is also present in the invariant mass spectrum. It was verified on Monte Carlo 
simulations that this does not bias the extracted signal yield, because the invariant mass distribution of 
these candidates is wide enough to be accounted for by the background function. Moreover, at low 
Pt, the particle identification selection strongly suppresses this contribution. For D*+ mesons, the mass 
difference Am distribution is fitted using a function that consists of a Gaussian describing the signal and 
the term a y^Am — niji ■ QH^"'-"'yt) for the background Q, where ntji is the charged pion mass, and a and 
b are free parameters. For all three D meson species, the mean of the Gaussian is compatible with the 
PDG value 1:211 within errors, and its width is well reproduced in the simulation. The extracted D meson 
raw yields are reported in Table [T] 

4 D meson cross sections 
4.1 Corrections 

The production cross sections of prompt charmed mesons were calculated as (e.g. for D+): 
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(GeV/c) (GeV/c) p^ (GeV/c) 

Figure 3: Acceptance x efficiency for D*^, D+, and D*+ mesons, as a function of pt (see text for details). 



raw^^^-j j-j^g measured inclusive raw yield, obtained from the invariant mass analysis in each pt 
interval (of width Apt); /prompt is the prompt fraction of the raw yield; (Acc x £)prompt is the acceptance 
times efficiency of prompt mesons, where e accounts for vertex reconstruction, track reconstruction and 
selection, and for D meson candidate selection with the secondary vertex and particle identification cuts 
described in section[3l Ay (= 2yfid) is the width of the fiducial rapidity coverage (see section I3l2l ) and BR 
is the decay branching ratio |[2TI . The factor 1/2 accounts for the fact that the measured yields include 
particles and anti-particles while the cross sections are given for particles only. The integrated luminosity 
was computed as Lint = App Ms/Cpp MB. where A'^pp MB and app.MB are the number and the cross section, 
respectively, of pp collisions passing the minimum-bias trigger condition defined in section|2l The CJpp mb 
value, 62.5 mb, was derived from a measurement using a van der Meer scan [25] of the cross section 
of collisions that give signals in both sides of the VZERO scintillator detector (cJvzero-and) ifTSl - The 
normalization factor, CJpp vzero-and / Cpp.MB ~ 0.87, was found to be stable within 1% in the data sample. 
The uncertainty on app.MB is determined by the systematic uncertainty of 3.5% on C7vzero-and> which 
is due to the uncertainties on the beam intensities ll26l and on the analysis procedure related to the van 
der Meer scan of the signal. 

The rapidity acceptance connection, using the factor ly^, with y^ varying from 0.5 at low pi to 0.8 
at high pi, assumes that the rapidity distribution of D mesons is uniform in the range \y\ < yM- This 
assumption was checked using the PYTHIA 6.4.21 event generator |[T9l with Perugia-0 tuning ll20l and 
the FONLL pQCD calculation ifTllTTl. both of which generate a D meson yield that is uniform within 1% 
in the range \y\ < 0.8. The (Acc x e) correction was determined using Monte Carlo simulations based on 
the GEANT3 transport code |[24l . The luminous region distribution and the conditions of all the ALICE 
detectors in terms of active channels, gain, noise level, and alignment, and their evolution with time 
during the 2010 LHC run, were included in the simulations. Proton-proton collisions were simulated 
using the PYTHIA 6.4.21 event generator |[T9l with Perugia-0 tuning 1201 . Only events containing 
D mesons were transported through the apparatus and reconstnicted, and the efficiency was extracted 
separately for prompt D mesons and D mesons from B meson decays. Figure [3] shows, as a function 
of transverse momentum, the acceptance times efficiency (Acc x e) for D*', D+, and D*+ mesons with 
bl < Jfid- At low Pt, the efficiencies are of order 1% or less, while for large pt the efficiencies increase and 
flatten at about 10-20% for D° and D+, and 30% for D*+. The efficiencies without particle identification 
selection, shown for comparison, are the same as those with particle identification, indicating that this 
selection is essentially fully efficient for the signal. The efficiencies for D'' and D+ mesons from B meson 
decays, also shown for comparison, are larger by about a factor of two. This behaviour is due to the fact 
that feed-down D mesons decay further from the primary vertex, because of the large B meson lifetime 
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(cT ^ 500 /im [211). For D*+ mesons, the efficiency for the prompt and feed-down components are the 
same in the pi range above 4 GeV / c, where no strong cuts on the separation of the D" decay vertex from 
the primary vertex are applied. 

The fraction /prompt of D mesons coming from c quark hadronization, i.e. the coiTcction factor that 
accounts for the feed-down from B meson decays, was evaluated using the B production cross section 
from the FONLL pQCD calculation dllZl. which describes well beauty production at Tevatron |[3l and 
at the LHC ||5j|6l, and the B — > D decay kinematics from the EvtGen package [21] . The computed cross 
section for the feed-down component for each of the three D meson species was used, together with the 
Monte Carlo acceptance times efficiency (Acc x £)feed-down for D mesons from B decays (see Fig.O, to 
compute the expected feed-down contribution in the measured raw yields: 

/prompt = 1 - (A^°^ ""/A^°" (2) 

• ^y^p^ ■ (Acc X £)feed-down " BR • Lim . (3) 



with: 



from B raw 



J— -D^ from B 

2 ""fonll 



|>i<vfld dpt 



lvko.5 



The symbol of the ;?t-dependence (pi) is omitted in the formulas, for brevity. The resulting prompt 
fraction /prompt is shown in Fig. |4] (left-hand panel) by the soUd horizontal lines, for the case of D° 
mesons. The prompt fraction ranges between 80% and 90%, depending on the pi interval, these values 
being determined also by the different efficiencies for prompt and feed-down D mesons. In order to 
estimate the systematic uncertainty, the perturbative uncertainty on the FONLL beauty production cross 
section was considered, as well as an alternative way of using the FONLL calculation. The former 
contribution was obtained by varying the b quark mass and the factorization and renormalization scales 
as suggested in, e.g., ll29l . The alternative method consisted of computing the prompt fraction using the 
FONLL cross sections for prompt and feed-down D mesons (with B — )• D via EvtGen [i28il for the latter) 
and their respective Monte Carlo efficiencies: 



/t 



/ 

(Acc X £) feed-down 



prompt 



-D+ from B 



dp, 



1 + - 



\ 

|y|<0.5 



(Acc X £)prompt da^^, 



\ 



dp, 



(4) 



|y|<0.5 / 



The resulting prompt fraction is shown by the dashed horizontal lines in Fig.|4](left-hand panel). The full 
envelope of the uncertainty bands from the two methods, which is shown by the boxes in the figure, was 
taken as a systematic uncertainty. The uncertainty related to the B decay kinematics was disregarded, 
after verifying that the difference resulting from using the PYTHIA lT9l decayer instead of EvtGen |l28l| 
is negligible with respect to the FONLL B meson cross section uncertainty. 

The prompt fraction of D" mesons in the reconstructed yield was also estimated with a data-driven 
method based on the measured impact parameter distribution of D° meson candidates in each pi interval, 
as done previously by the CDF Collaboration 171. This method exploits the different shapes of the distri- 
butions of the impact parameter to the primary vertex of prompt and feed-down (displaced) D mesons. 
The impact parameter distribution of D" mesons was obtained from the one measured for candidates 
with invariant mass in the range — Mqo| < 2 a, after subtracting the background contribution esti- 
mated from the candidates in the side-bands (in the range 4.5 a < \m — Mpoj < 4.5 a + 100 MeV/c^). 
The prompt fraction was estimated by fitting the resulting impact parameter distribution with a two- 
component function. The first component is a detector resolution term, modelled by a Gaussian and 
an exponential term, describing the impact parameter of prompt D mesons. The second component ac- 
counts for the reconstructed impact parameter distribution of D from B decay, which is modelled by a 
convolution of the same detector resolution term with a double-exponential function describing the true 
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Figure 4: Left; prompt fraction /prompt of the D*' raw yield as a function of pu for the two FONLL-based methods 
(solid: central value, from Eq. (|2]i; dashed: alternative method, from Eq. Q) and for the impact parameter fit 
method (circles); the boxes show the envelope of the uncertainty bands of the two FONLL-based methods; the error 
bars show the total uncertainty from the impact parameter fit, including the statistical and systematic contributions. 
Right: an example of meson impact parameter distribution in the transverse plane, for 2 < pt < 3 Ge V / c; 
the distribution is background-subtracted and fitted with the two-component function for prompt and feed-down 
contributions, as described in the text; the resulting prompt fraction, impact parameter resolution for prompt 
mesons, and ;i;'/(number of degrees of freedom) of the fit are given; the inset, with linear scale, shows also the 
negative entries, resulting from the background subtraction. 



impact parameter of secondary D mesons. The fit parameters are the width of the Gaussian and the frac- 
tion of prompt D mesons, that is, the relative weight of the prompt and secondary D meson components. 
An example of such a fit is shown in Fig. |4] (right-hand panel) for the D° mesons in the transverse mo- 
mentum interval 2 < pi <3 GeV/c. The prompt fraction of D° mesons, measured with this method for 
2 < Pt< 12 GeV/c, is shown by the circles in Fig.|4](left-hand panel) and is found to be in general agree- 
ment with the FONLL-based estimations. Because of the large background for pt <2 GeV/c and the 
poor statistics available for pt> 12 GeV/c, this method was used only as a check of the FONLL-based 
prompt fraction estimation, for the pt intervals with large signal yield. 

4.2 Systematic uncertainties 

Several sources of systematic uncertainty were considered, including those affecting the signal extraction 
from the invariant mass spectra and all the correction factors applied to obtain the ;?t -differential cross 
sections. A summary of the estimated relative systematic uncertainties is given in Table [2l for the lowest 
and highest pt interval (see Table [U for each meson species. All the uncertainties evolve monotonously 
as a function of pt, except for the yield extraction and feed-down correction contributions, which ai^e 
larger at low and high pi than at intermediate pi. 

The systematic uncertainty on the yield extraction from the invariant mass spectrum in a given pi interval 
was determined by repeating the fit in a different mass range, and also varying the function to describe 
the background. A polynomial, instead of an exponential, was used for D'' and D+ mesons, while a 
power law convoluted with an exponential and a polynomial was considered for D*+ mesons, instead of 
the function defined in section ll!2l A method based on bin counting (after subtraction of the background 
estimated from a fit in the mass side bands) was also used. The uncertainty was defined as the maximum 
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variation in the extracted yields from these different methods. 

The systematic uncertainty related to the tracking efficiency includes the effects arising from track finding 
in the TPC, from track propagation from the TPC to the ITS, and from track quality selection. It was 
estimated from the comparison of data and simulation and from the variation of the track selection. The 
resulting uncertainty is 8% for the two-body decay of D*' mesons and 12% for the three-body decay of 
D+ mesons. For the D*+ case, a slightly larger systematic uncertainty of 13% was assigned in the low 
Pi region (below 3 GeV/c), because the soft pion often has pi< 150 MeV/c and is reconstructed only in 
the ITS. The tracking efficiency in this case has a significant uncertainty arising from the description of 
hadronic interactions in the simulation of the detector response. 

A systematic effect can arise due to residual discrepancies between data and simulation for the variables 
used to select the signal D meson candidates. The distributions of these variables were compared for 
candidates passing loose topological cuts, i.e. essentially for background candidates, and found to 
be well described in the simulation. The systematic effects due to residual differences between data 
and simulation were quantified by repeating the analysis with different sets of cuts. In particular, the 
cut values were changed in order to vary the signal by at least 20% below pt = 8 GeV/c. From the 
corresponding variation of the corrected spectra, a systematic uncertainty of about 10% was estimated 
for each D meson species. As a further cross-check, the secondary vertices in the simulation were 
reconstructed also after a track-by-track scaling by a factor 1.08 of the impact parameter residuals 
with respect to their true value. This scaling is aimed at reproducing the impact parameter resolution 
observed in the data (see Fig. [T]) and accounts for possible residual detector misalignment effects not 
fully described in the simulation. The resulting variation of the efficiency was found to decrease from 
4% at pi = 1-2 GeV/c to less than \% for p^ > 5 GeV/c. This effect was not included explicitly in 
the uncertainty estimation, since it is to some extent accounted for in the cut variation study and its 
magnitude is much smaller than the 10% uncertainty assigned to the cut efficiency corrections. 

The systematic uncertainty induced by a different efficiency for particle identification in data and 
simulation was evaluated by repeating the analysis either without applying this selection, or with a tighter 
selection (2 a compatibility instead of 3 a). The variation of the corrected yields obtained without PID 
and with these two selections (standard and tighter) was assigned as a systematic uncertainty. Typical 
values are 3-5%, with the exception of 15% for D+ in 1-2 GeV/c. 

The accuracy of the description of the evolution of the experimental conditions with time was verified 
by analyzing separately sub-samples of data collected during different periods and with different orien- 
tations of the magnetic field. The results were found to be compatible within statistical uncertainties for 
all three meson species. Furthermore, the /^t-differential yields for each D meson measured separately 
for particles and anti-particles were found to be in agreement within statistical uncertainties. 

The effect of the shape of the simulated D mesons spectrum within our pt intervals was estimated from 
Table 2: Summary of relative systematic uncertainties for the lowest and highest interval for each meson species. 
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the relative difference in the Monte Carlo efficiencies obtained with the pt shapes from PYTHIA |[T9l 
with Perugia-0 tune |[20l and from the FONLL pQCD calculation |[T][33. These two models predict a 
significantly different slope at high pt (dN /dpt oc p:^'^-^ for FONLL and <x p^^-^ for PYTHIA Perugia- 
0), which however results in a systematic effect on the D meson selection efficiency of only 3% for 
I < Pi <2 GeV/c, and less than 1% at higher pi. 

The systematic uncertainty from the subtraction of feed-down D mesons from B decays was estimated 
as described in the previous section. It ranges between +45% at low pt (1-2 GeV/c) and +10% at high pt 
12 GeV/c). 

Finally, the results have global systematic uncertainties caused by the branching ratios (taken from ^TDi ) 
and by the minimum-bias cross section (3.5%). 



4.3 Results 



The /Jt-differential inclusive cross sections ai^e shown in Fig. [5] for prompt D'', D+, and D*+ mesons. 
The error bars represent the statistical uncertainties, while the systematic uncertainties are shown as 
boxes around the data points. The numerical values for the differential cross sections are reported in 
Table [3l together with their statistical and systematic uncertainties, as well as the average transverse 
momentum (pi) of D mesons in each pi interval. The value of {pi) was obtained from the meson pi 
distribution in the considered interval, after subtracting the background contribution from the side bands 
in the invariant mass distribution. The resulting (pt) values for the three D meson species are compatible 
within uncertainties in all pi bins and their average is reported in the table. 

The measured D meson inclusive differential production cross sections are compared to two theoretical 
predictions, namely FONLL 1111301 and GM-VFNS El HI. Both calculations use CTEQ6.6 parton 
distribution functions (PDF) Il32l and vary the factorization and renormalization scales, /ip and /Ir, 
independently in the ranges 0.5 < /ip/wit < 2, 0.5 < jJ-R/nii < 2, with the constraint 0.5 < /Xf/Mr < 2, 
where nit = \/pf+ni^- The chaim quark mass is varied in FONLL within 1.3 < rric < 1.7 GeV/c^. 
Both calculations are compatible with the measurements, within the uncertainties. The central value of 
the GM-VFNS predictions lies systematically above the data, while that of the FONLL predictions lies 
below the data. For FONLL, this feature was observed also at ^/s = 0.2 TeV (pp) lllOllllll and 1.96 TeV 
(pp) lEl- With a reach down to /^t = 1 GeV/c, this measurement probes the gluon distribution in the 
X range of a few 10^^. Within the current uncertainties of the experimental measurement and of the 



Table 3: Production cross section in |y| < 0.5 for prompt D'', D+, and D*+ mesons in pp collisions at ^/s = 7 TeV, 
in transverse momentum intervals. The normalization systematic uncertainty (3.5% from the minimum-bias cross 
section plus the branching ratio uncertainties, as of Table|2]i is not included in the systematic uncertainties reported 
in the table. 
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Figure 5: (colour online) pt-differential inclusive cross section for prompt D", D+, and D*+ mesons in pp 
collisions at = 7 TeV compared with FONLL |[l]|30l and GM-VFNS gUni theoretical predictions. The 
symbols are positioned horizontally at the centre of each pi interval. The normalization uncertainty is not shown 
(3.5% from the minimum-bias cross section plus the branching ratio uncertainties, as of Table|2]). 



theoretical predictions, it is not possible to draw conclusions about small-x gluon saturation effects (see 
section [T]). 

The /^t-integrated visible cross sections, o^^^{pi > 1 GeV/c, \y\ < 0.5), for the three mesons were 
extrapolated down to /^t = to estimate the production cross sections per unit of rapidity da/dy at 
mid-rapidity. The extrapolation factor was computed from the FONLL calculation |[T]|221 as the ratio 
(daFONLL/dj)/a^gNLL and it amounts to 1.25+°;^^ for D° and D+, and 1.21+°;29 for D*+, for the central 
values of the calculation. Its uncertainty was obtained as a quadratic sum of the uncertainties from 
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charm mass and perturbative scales, varied within the aforementioned ranged, and from the CTEQ6.6 
PDF sets [|32J . The cross sections for the three mesons are: 

da^Vdj = 516±41(stat.)li75(syst.) ± 18(lumi.) ±7(BR)l'^°(extr.) jxh, 

da°^/d3; = 248±30(stat.)+92(syst.)±9(lumi.)±5(BR)+^7(extr.) ^h, 
da^'^/dy = 247 ±27(stat.)^8i(syst.) ±9(lumi.) ±4(BR)l^^(extr.) jUb. 

5 Summary 

We have presented the measurement by the ALICE Collaboration of the inclusive differential production 
cross sections of prompt D mesons at central rapidity, in pp collisions at ^/s = 7 TeV within I < pt < 
24 GeV/c. D mesons were reconstructed in the decay channels D" — ^ K^;r+, D+ — )• K^n^n^, and 
D*+ — D°;r+, and their charge conjugates. 

The /7t-differential cross sections are reproduced within uncertainties by theoretical predictions based on 
perturbative QCD, FONLL 1301 and GM-VFNS EB- More in detail, the data tend to be higher than the 
central value of the FONLL predictions, as it was observed also at lower collision energies, at RHIC and 
at the Tevatron |l7l[T0l[TTl. For GM-VFNS, instead, the data lie on the lower side of the predictions, at 
variance with the case of Tevatron energy QIH, indicating that the energy dependence is steeper in this 
model than in data. Our results, together with existing measurements at lower energies, can contribute to 
a better understanding of charm production in pQCD. 

Furthermore, the measurements that we have presented open new possibilities to test PDF dynamics, 
in the regime of paiton fractional momentum below x ~ 10^"^ and squared momentum transfer down 
to ~ (4 GeV)^, where the onset of gluon PDF saturation effects has been conjectured |[T2l . Within 
the uncertainties of the data and of the theoretical predictions, the framework of factorized QCD cal- 
culations provides a reasonable description of the data points down to the lowest measured transverse 
momentum. However, accurate calculations incorporating saturation effects are needed in order to draw 
firm conclusions on their relevance for low-momentum chaim production at LHC energies. 
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